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Abstract 

A dynamic scaling Ansatz for the approach to the Self-Organized Critical (SOC) regime is 
proposed and tested by means of extensive simulations applied to the Bak-Sneppen model (BS), 
which exhibits robust SOC behavior. Considering the short-time scaling behavior of the density of 
sites (p(t)) below the critical value, it is shown that i) starting the dynamics with configurations 
such that pit = 0) — > one observes an initial increase of the density with exponent 9 = 0.12(2); 
ii) using initial configurations with pit = 0) — * 1, the density decays with exponent 5 = 0.47(2). It 
is also shown that he temporal autocorrelation decays with exponent C a = 0.35(2). Using these, 
dynamically determined, critical exponents and suitable scaling relationships, all known exponents 
of the BS model can be obtained, e.g. the dynamical exponent z = 2.10(5), the mass dimension 
exponent D = 2.42(5), and the exponent of all returns of the activity tall = 0.39(2), in excellent 
agreement with values already accepted and obtained within the SOC regime. 
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Nearly two decades ago, Bak et al. .1] proposed the celebrated concept of Self-Organized 
Criticality (SOC) in order to describe complex systems capable of evolving toward a critical 
state without the need of tuning any control parameter. This is in contrast to the case 
of standard critical behavior where critical points are reached by tuning a suitable control 
parameter (temperature, pressure, etc.). The study of SOC behavior has attracted huge 
attention due to its ubiquity in a great variety of systems in the fields of biology (evolutionary 
models), geology (earthquakes), physics (flick noise), zoology (prey-predators and herds), 
chemistry (chemical reactions), social sciences (collective behavior of individuals), ecology 
(forest-fire), neurology (neural networks), etc. 0,0]. 

In spite of the considerable effort invested in the study of SOC, it is surprising that little 
attention has been drawn to the understanding of the dynamic approach to the SOC regime 
when a system starts far from it. This issue is relevant for a comprehensive description 
of the phenomena since the SOC state behaves as an attractor of the dynamics. For the 
case of standard criticality, the existence of a short-time universal dynamic scaling form 
for model A has only recently been established. In fact, according to a field-theoretical 
analysis followed by an e— expansion m, which was subsequently been extensively confirmed 
by means of numerical simulations p, a short-time universal dynamic evolution that sets 
in right after a time scale t m i C , which is large enough in the microscopic sense but still very 
small in the macroscopic one, has been identified. It is worth mentioning that by means 
of short-time measurement one can not only evaluate the dynamic exponent z and relevant 
(static) exponents, but also the exponent (0) describing the scaling behavior of the initial 
increase of the order parameter of model A. 

Within this context, the aim of this work is to propose a dynamic scaling Ansatz to 
describe the approach to the SOC state. Furthermore, our proposal is validated by extensive 
simulations of the evolutionary Bak-Sneppen model, showing that the dynamic approach to 
the SOC state is in fact critical and its study allows us to evaluate exponents that are in 
excellent agreement with independent measurements performed within the SOC regime. 

The Bak-Sneppen (BS) model is aimed at simulating the evolution of life through individ- 
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ual mutations and their relation in the food chain jy, LD, |8|, |9( . Each site of a d— dimensional 
array of side L represents a species whose fitness is given by a random number / taken from 
a uniform distribution P(f) in the range [0, 1]. The system evolves according to the follow- 
ing rules: (1) The site with the smallest fitness is chosen. (2) A new fitness is assigned to 
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that site, i.e., a random number taken from P. This rule is based on the Darwinian survival 
principle, i.e., the species with less fitness are replaced or mutated. (3) At the same time, 
the fitness of the nearest-neighbor sites are changed. This rule simulates the impact of the 
mutation over the environment. 

The BS model is the archetype example of extremal dynamics and perhaps, it is the 
simplest example of a system exhibiting a robust SOC behavior. We will perform simulations 
in d — 1, assuming periodic boundary conditions. It is well known that the system reaches 
a stationary (SOC) state where the density p of sites with fitness below a critical value is 
negligible (/ < f c , with f c = 0.66702(3)), but it is uniform above f c p. Within the SOC 
state, the BS model exhibits scale-free evolutionary avalanches and punctuated equilibrium 



The dynamic scaling Ansatz. The short-time dynamic scaling form has originally been 
formulated for the Ising model with two state spins p. In the absence of magnetic fields, 
the Ising magnet exhibits a second-order phase transition when the temperature is tuned 
at the critical point (i.e., standard critical behavior). By analogy, we also considered that 
in the BS model, each site can only be in one of two possible states af. occupied (erj = 1) 
when its fitness is below f c or empty ((Tj = 0) when its fitness is above f c . Furthermore, 
since the magnetization goes to zero at the critical temperature of the Ising system, while in 
the BS model the density of sites p also vanishes when the SOC regime is reached [jj, it is 
also reasonable to propose a scaling Ansatz for the density. However, it is worth mentioning 
that the BS model lacks a control parameter, such as the temperature for the Ising model. 
Hence, in the limit L — > oo, the proposed scaling reads 



where po is the initial density, z is the dynamic exponent, xq is the exponent of the rescaling 
of po> N is an exponent, and b is a scaling variable. 

Simulation results and discussion. Computer simulations were performed in ensembles 
of 10 3 different systems having the same initial density of sites with fitness below f c . Notice 
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FIG. 1: Temporal evolution of average density for several values of initial density pq as listed in 
the figure. 

that usually the site with the smallest fitness value at time t is called the active site. Figure ^ 
shows the temporal evolution of density p(t, p ) as obtained for different values of p . Three 
different regimes can be distinguished: i) A short-time regime (0 < t < t rnax « 10 3 ), 
which holds for low initial densities pn < 0.1, where the density exhibits an initial increase. 

n 

Since the value of the "effective" exponent slightly depends on p , as ussually p, we have 
performed an extrapolation obtaining 9 = 0.12(2) in the po — ► limit, ii) An intermediate- 
time regime (t max < t < t stat « 10 6 ) where the density decreases also following a power-law 
with exponent 5 = 0.47(2), and iii) A long-time regime (tstat < t)-> where the system arrives 
at a stationary state with a constant average density p s tat- 

In order to understand the observed behavior it is useful to analyze first a simulation 
started with a single site (po = l/£)- One observes (not shown here for the sake of space) a 
single avalanche and the density increases monotonically with exponent 9, as in the case of 
figure ^ for p Q —>■ 0. The spatiotemporal evolution of the avalanche is delimited according to 



where D = 2.43(1) |2| is the mass dimension exponent. 

Then it is convenient to consider the evolution of n e epidemics started simultaneously 
with a separation of r e empty sites between them. Two neighbor epidemics collide when 
each of them expands its activity over y sites, on average. Then, the border between them 
disappears leading to a single epidemic. So, according to equation (J2J), this requires a time 



r(t) ~ t 



l/D 



(2) 
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of order t e ~ (r e /2) D , hence for the collision of n e epidemics the time required is of order 
t ~ n e t e . Then, one expects an initial increase of the density until t max , given by 

tmax ~ n e j . (3) 

Now, for an initial random distribution, one has that n e = Lp and the distance between 
particles at t = is of order r e ~ 1/po- Then, equation ® becomes 

tmax ~ PI~ D - (4) 

1—1 

Also, the system reaches the stationary state for a time of order t s t a t ~ L D [9]. In the 
thermodynamic limit p stat — > and t sta t oo. 

It is worth mentioning that the time t used in this work is a discrete sequential time. 
An alternative definition corresponds to the parallel time t», usually employed to define the 
dynamic exponent z according to 

*||~r«. (5) 

The unit of parallel time is defined as the average number of actualization steps that have 
to be performed to change the state of all the occupied sites (n(t)), then 

t-*t+l and t\\ -> tii + -4t (6) 

n(t) 

From these definitions one has tn ~ t 1 ' 6 , which inserted in equations (0) and (J2J), yields 

The time scaling behavior of the density can be obtained from equation (P) by replacing 
b = t 1 ^, which yields 

p(t,p Q )=t-«/*$(t*°/*p ) (8) 

where $ is a scaling function. For po ~~ * 1 an d within the intermediate time regime, p{t) 
becomes independent of p (see figure HJ), thus assuming x /z > one has $(t X0//z p ^> 1) = 
constant, which holds for t ^> t max oc Pq Z ^ x °. Then, t max sets the time scale for the initial 
increase of the density, as shown in figure ^ [nj . Furthermore, according to equation (j3J) 
the following relationship between exponents should hold 

D — 1 = z/x . (9) 
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FIG. 2: (a) Log-log plots of p versus po obtained for times t = 10,20,50, 100,200, and 500 (from 
bottom to top), (b) Scaling plot of the curves shown in (a) using equation (|lljl. 



Also, p should decrease according to 



p(t) ~ t- 



(10) 



with 5 = K/z = 0.47(2). On the other hand, for p — > and within the short-time regime 
(t <C t max ,t — > 0) the dependence of $ on the initial density becomes relevant, thus we 
assume $(x) ~ x u , where u is an exponent. Hence, replacing in equation (JBJ), one has 



p(t,p ) =plt\ 



(11) 



where 9 = —5 + u—. So, = 0.12(2) is the exponent that describes the initial increase of 
the density within the short-time dynamic critical regime of the BS model. 

In order to calculate u, the dependence of the density on p (i.e., p(po)) was measured 
at different times, as shown in figure El (a). Also, the scaling Ansatz suggested by equation 
(llljl is shown in figure El (b). The observed data collapse is satisfactory and the exponent 
u = 0.87(3) was measured. 

On the other hand, replacing 21 = in equation (|BJ) one obtains 



p(t,p ) = r 8 *(j> t 9 ?). (12) 

The shape of the scaling function $ is shown in figure El and the excellent data collapse 
obtained by plotting the data already shown in figure ^ strongly supports the formulated 
scaling hypothesis and the calculated exponents. In fact, the best fit of the data obtained, 
for p - * 0, gives u = 0.89(3) (see dotted line covering two decades in figure EJ) that is in 
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agreement with the preliminary estimation already performed with the data shown in figure 
121 (covering less than one decade). 

Also, the scaling of the initial density can also be obtained by replacing b = p 1 ' Xo in 
equation (JTJ), giving 

p(t,po)=Po /x °n(tpo /XO ), (13) 

where Q is a scaling function that, for t — > and po ~~ 0, can be approximated by Q(x) — > 
(a; — ► 0), where t> is an exponent. Thus one obtains 

V,/xo+v — 

p(t,Po)= Po x °t v . (14) 

Now, by comparing equation dHj) with equation (fTTjl it follows that v = 9 and u = ^ + 
Hence, 

^ = ^ = (15) 

Xo Xo/-2 5 + 

Then, using ^ = ^ and equation (|THj) . one has that equation (JT3J) can be written in 
terms of already measured exponents, so 

uS u 

p(t, Po )= P r^(tpr)- (is) 

Figure 0] shows plots of the numerical data performed according to equation (|16)l . The 
shape of the scaling function Q can be observed and the collapse of the curves also supports 
the formulated scaling hypothesis. 

Valuable information on the dynamic behavior of the system can also be obtained by 
measuring the temporal autocorrelation of the state of the site (A(t,to), averaged over all 
sites) that is expected to decay according to a power law jj, 0|, namely 

A(t Q ,t) = {(Ti(jt )(7i(t)) - (<7 4 (* )) <<*(*)) ~ t~° a (17) 

where cjj(t) = 1,0 is the state of the site at time t. Also [^| 

C a = d/z - 6. (18) 
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FIG. 3: Scaling plot of the data shown in figure ^ obtained by using equation (|12j) for the time- 
scaling behavior. The dotted line -slightly shifted up for the sake of clarity- with slope u = 0.89 
corresponds to the best fit of the data obtained for the lowest density. 





FIG. 4: Scaling plot of the data shown in figure ^ obtained by using equation (|16|) for the scaling 
of the initial density. 

Figure El shows log-log plots of the autocorrelation versus t obtained using very low initial 
densities since equation (|T7|) is expected to hold for p n —* 0. From these plots one determines 
C a = 0.35(2), and replacing this value in equation (|T5|) the dynamic exponent z = 2.13(5) 
is obtained. Then by using equation (J7J) we obtain D = 2.42(5). These exponent values are 
in agreement with those already published in the literature that were obtained within the 
SOC regime jjj, i.e., z = 2.1(5) and D = 2.43(1). Furthermore, equation (jHJ) can be tested, 
yielding D — 1 = 1.42(5) (left-hand side) and z/xq — u/{9 + 5) — 1.47(5) (right-hand side). 

Now, we can discuss the relationship between short-time dynamic measurements and the 
behavior of relevant distributions characterizing return times of the activity to a given point 
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in space. In fact, the distribution Pall(Q, t\\) is the probability for the activity at time tn 
to revisit a site that was visited at time and is often referred to as the distribution of 
all return times jf) . For tn 3> 1 such a distribution decays as a power law of the form 
Pall{^i t^ TALL , where ta,ll is the 'lifetime' exponent for all returns of the activity |9j. 
Then, considering that A(0, t) is the probability that a given site being occupied at t = 
returns to be occupied at time t, one has that A(Q,t) oc Pall(0, ty) since a site can change 
its state only when such a site, or any of its two neighbours, are visited by the activity. 
Then, by recalling that we are working in terms of the discrete sequential time t such as 
t\\ oc t l ~ e Q, one has 

A(t, 0) ~ t-d-^", with tall = (19) 

1 — u 

By means of short-time dynamic measurements we have obtained tall — 0.39(4), in 
excellent agreement with the value calculated within the stationary state (t > 10 6 ), given 
by r ALL = 0.42(2) 0. 

Summing up, the proposed dynamic scaling Ansatz for the approach to the stationary 
state generalizes the concepts previously developed to describe the critical dynamics of model 
A for systems exhibiting SOC. The dynamic scaling behavior was tested with the BS model 
showing that it holds for the density of sites with fitness below the critical one. Remarkably 
the exponents calculated by using the dynamic approach are in excellent agreement with 
those already measured within the stationary state. Using well established relationships 
between exponents and considering that only two of them are basic, the dynamic measure- 
ment allows the self-consistent evaluation of all exponents of the BS model. So, we have 
shown that the dynamic approach to the SOC regime is indeed critical and it is governed 
by the dynamic exponent z, standard exponents (i.e., those available from stationary deter- 
minations), and by the exponent 9 that is introduced to describe the initial increase of the 
density. It is worth emphasizing that the criticality observed in the approach to the SOC 
state not only addresses new and challenging theoretical aspects of SOC behavior, but also 
is of great practical importance for the determination of relevant exponents. 
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FIG. 5: Log-log plots of the temporal autocorrelation versus t obtained for different values of pq. 
The dotted line with slope C a = 0.35 has been drawn for the sake of comparison. 
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